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Abstract 

The classical (i.e. non-quantum) equilibrium statistical mechanics of a two-dimensional 
one-component plasma (a system of charged point-particles embedded in a neutralizing back- 
ground) living on a pseudosphere (an infinite surface of constant negative curvature) is con- 
sidered. In the case of a flat space, it is known that, for a one-component plasma, there are 
several reasonable definitions of the pressure, and that some of them are not equivalent to 
each other. In the present paper, this problem is revisited in the case of a pseudosphere. 
General relations between the different pressures are given. At one special temperature, the 
model is exactly solvable in the grand canonical ensemble. The grand potential and the 
one-body density are calculated in a disk, and the thermodynamic limit is investigated. The 
general relations between the different pressures are checked on the solvable model. 
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1 INTRODUCTION 



Coulomb systems such as plasmas or electrolytes are made of charged particles interacting 
through Coulomb's law. The simplest model of a Coulomb system is the one-component plasma 
(OCP), also called jellium: an assembly of identical point charges, embedded in a neutralizing 
uniform background of the opposite sign. Here we consider the classical (i.e. non-quantum) equi- 
librium statistical mechanics of the OCP. Although many features of more realistic systems are 
correctly reproduced, this model has the peculiarity that there are several reasonable definitions 
of its pressure, and some of these definitions arc not equivalent to each other [1, 2]. 

The two-dimensional version (2D OCP) of the OCP has been much studied. Provided that 
the Coulomb potential due to a point-charge is defined as the solution of the Poisson equa- 
tion in a two-dimensional world (i.e. is a logarithmic function — Inr of the distance r to that 
p oint- charge ) , the 2D OCP mimicks many generic properties of the three-dimensional Coulomb 
systems. Of course, this toy logarithmic model does not describe real charged particles, such 
as electrons, confined on a surface, which nevertheless interact through the three dimensional 
Coulomb potential 1/r. One motivation for studying the 2D OCP is that its equilibrium sta- 
tistical mechanics is exactly solvable at one special temperature: both the thermodynamical 
quantities and the correlation functions are available [3]. 

How the properties of a system are affected by the curvature of the space in which the 
system lives is a question which arises in general relativity. This is an incentive for studying 
simple models. Thus, the problem of a 2D OCP on a pscudosphere has been considered [4]. A 
pseudosphere is a non-compact Riemannian surface of constant negative curvature. Unlike the 
sphere it has an infinite area and it is not embeddable in the three dimensional Euclidean space. 
The property of having an infinite area makes it interesting from the point of view of Statistical 
Physics because one can take the thermodynamic limit on it. 

For this 2D OCP on a pseudosphere, the problem of studying and comparing the different 
possible definitions of the pressure also arises. This is the subject of the present paper. In 
Section 2, we give some basic properties of the pseudosphere and of a 2D OCP on it. In Section 3, 
we define the different pressures and derive general relations between them. In Section 4, we 
illustrate the general properties by considering the special temperature at which all properties 
can be explicitly and exactly calculated. 



2 PSEUDOSPHERE AND ONE-COMPONENT PLASMA 
2.1 The Pseudosphere 

There are at least three commonly known sets of coordinates to describe a pseudosphere of 
Gaussian curvature — The one which renders explicit the resemblance with the sphere is 

Q = (<Z^) Q^) = (^^) Q"^) = (''") ^) with r G [0, oo[ and (p G [0, 2tt[, the metric being 

ds^ = ga(3 = a^{dT^ + sinh^ r dip"^) (2.1) 

Another set of coordinates often used is {r,ip) with r/(2a) = tanh(r/2). They are the polar 
coordinates of a disk of radius 2a. The metric in terms of these new coordinates is 

- [1 - (rV4a2)]2 ^^-^^ 

The disk with such a metric is called the Poincare disk. Through an appropriate conformal 
transformation, the Poincare disk can be mapped onto the upper half-plane with some metric, 
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the Poincare half-plane, but this latter representation will not be used here. The geodesic 
distance doi between any two points % = (tq, (po) and qi = (ri, (pi) on the pseudosphere is given 

by 

cosh((ioi /o) = cosh ti cosh tq — sinh ri sinh tq cos((^i — (^o) (2-3) 

Given the set of points at a geodesic distance from the origin less than or equal to d, that 
we shall call a disk of radius d, we can easily determine its circumference 

C = 27r o sinh {-\ /y^ vr a e*^/" (2.4) 



and its area 



^ = 47ra^sinh2(£) (2-5) 



The Laplace-Beltrami operator on the pseudosphere is 

A- — f^^— sinhr— (2 6) 

\ sinh r dr dr sinh^ r dip"^ J 

2.2 The One-Component Plcisma 

The 2D OCP which is considered here is an ensemble of N identical point particles of charge q, 

constrained to move in a disk of radius d = utq by an infinite potential barrier on the boundary 
of this domain r = tq. The average particle number density is ra = N/A, where A is the area 
(2.5). There is a background with a charge density pb = —qn^ uniformly smeared on the disk 
{pb is outside the disk). It is convenient to introduce the number of elementary charges in the 
background: Nb = UbA. The total charge is not necessarily 0, thus in general Ub 7^ n. 

The pair Coulomb potential v{d) between two unit charges, a geodesic distance d apart, 
satisfies the Poisson equation on the pseudosphere, 

Ai;(c/) = -27r(5(2)((i) (2.7) 

where (5^^^ {d) is the Dirac delta function on the curved manifold. This Poisson equation admits 
a solution vanishing at infinity. 



v{d) = — In 



tanh — 



(2.8) 



2a, 

The electrostatic potential of the background w{q) satisfies 

Au;(g) = -2'Kpb (2.9) 

By symmetry, this electrostatic potential is only a function of r. Expressing the Laplacian (2.6) 

in terms of the variable coshr, and requesting the solution to be regular at r = and to have the 
correct value at r = tq (corresponding to the background total charge), one finds the solution 



w{t) = 2TTa^qnb { In 



1 -tanh2(ro/2) 



+ sinh2(To/2) ln[tanh2(ro/2)] \ (2.10) 



1 -tanh2(r/2) 

Let dS = 27ra^ sinh rdr be an area element. The self energy of the background is 

= \ j pbw{T)dS 

^ J r<TO 

= (27ra2gn6)2{sinh2(ro/2) - ln[cosh2(To/2)] - sinh'^(ro/2) ln[tanh2(To/2)]} (2.11) 
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The total potential energy of the system is 

U = vo + Vpb + Vpp (2-12) 
where Vpp is the potential energy due to the interactions between the particles, 

1 ^ 

Vpp= 2l2^'^'"(^v^ (2.13) 
and Vpb is the potential energy due to the interaction between the particles and the background, 

N 
i=l 

3 THE DIFFERENT PRESSURES AND THEIR RELATIONS 

In the case of a flat system, the pressure which is often considered, termed the thermal pressure, 
is defined from the free energy F by the standard relation P^^^ = —{dF/dA)i3^N,Nt^ where P is 
the inverse temperature. In the case of a flat neutral (iV = Nf,) 2D OCP. this thermal pressure 
is given by the simple exact expression /3pW = n[l-{pcf/4:)] [5,6]. Thus, this thermal pressure 
becomes negative for /3g^ > 4, i.e. at low temperatures. This pathology of the OCP occurs also 
in three dimensions; it is related to the presence of an inert background without kinetic energy. 
Indeed, the uniform background can be considered as the limit of a gas of negative particles of 
charge — e and number density ly, when e ^ 0, ^ oo, for a fixed value of the charge density 
—ev. In this limit, the ideal-gas part (kinetic part) of the background average energy density 
becomes infinite. In the OCP Hamiltonian, this infinite energy density is omitted. The price 
paid for this omission is that the corresponding (infinite) ideal-gas contribution to the pressure 
is omitted, and the remaining pressure may be negative.^ 

Unhappy with this negativeness, Choquard et al.[l] and Navet et al.[2] have introduced 
another pressure, the kinetic pressure P^'^\ which is the pressure exerted on the wall by the 
particles of charge q only. This kinetic pressure turns out to be also the one which is obtained 
through the use of the virial theorem. Although for usual fluids the thermal and kinetic pressures 
are equivalent, in the presence of a background they are different, with the kinetic pressure being 
always positive. This positiveness led the above authors to argue that the kinetic pressure is the 
"right" one. Anyhow, a detailed comparison of the diverse possible definitions of the pressure 
of a flat OCP has been done [1]. 

In the present paper, it is this comparison that we extend to the case of a 2D OCP on a 
pseudosphere. We shall restrict ourselves to the case of a domain in the shape of a disk. We are 
especially interested in the thermodynamic limit, i.e. when the disk radius becomes infinite, for 
fixed values of /3, n, rz;,. 

3.1 Kinetic and Virial Pressures 

The average force exerted by the particles on a perimeter element ds is {1 / (3)n^^\To)ds , where 
n^^^ (r) is the one-body density at the distance ar from the origin. Therefore, the kinetic pressure 

''in the case of a two-dimensional iwo-conaponent plasma made of point-particles, the pressure also becomes 
negative when extrapolated to low temperatures fSq^ > 4. However, now /3q^ > 4 is outside the domain of 
definition of the partition function. 
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IS 

p('=) = (l//3)nW(To) (3.1) 

Wc assume that this quantity has a hmit when tq ^ oo. In Section 4, this assumption wih be 
checked in the special case (5q^ = 2. It will now be shown that the virial pressure P^""^ , i.e. the 
pressure computed from the virial theorem, is the same as P^^\ 

In terms of the 2N coordinate components and 2N conjugate momentum components 

, the Hamiltonian of our OCP of N particles is 

H{q^,p^)=T{q^,p^) + U{q'') (3.2) 
where U = U +confining potential and the kinetic energy T is 

1 ^ 

T=^T.9'"H^)PiaPif3 (3-3) 

i=l 

The Roman indices label the particles, and the lower or upper Greek indices denote covariant or 
contravariant components, respectively. As usual, a sum over repeated Greek indices is tacitly 
assumed. The equations of motion for particle i are 

BM 1 

(3.4) 

dH 1 dgf^^ dU 



- dqf ~ 2m dq^ ^'"^'^ Oq 



a 



where the dot stands for total derivative with respect to time. If we take the time derivative of 
Ei QiPiT = Ei nPiT, we find ^ 

-nZ^nPiT = — y,9 ^[qi)PiTPii3 - 7;—2^n^ — PipPi-f -}^n^ (3.5) 

at ^ m -Hf 2m ^ OTi -Hf oti 

t 1=1 1=1 1=1 

where the last term is called the virial of the system. Since the system is confined in a finite 
domain, the coordinates Ti{t) and their canonically conjugated momenta (0 remain finite at 
all times. Thus, 

{j,Y.nP.rJ=0 (3.6) 

where (. . .)t denotes a time average. Assuming that the system is ergodic, we can replace time 
averages by microcanonical averages. Assuming the equivalence of ensembles in the thermo- 
dynamic limit, wc can as well use canonical or grand-canonical averages (...). In the present 
section, we use canonical averages. The average of the r.h.s. of (3.5) vanishes. Separating in the 
last term of (3.5) the contribution from the forces exerted by the walls, which is, in the average, 
we obtain 



m . , 

1 1=1 



®One may be tempted to start with the time derivative of ^^qfpia = '^iinpiT + <fiiPiip). Note however 
that this quantity docs not remain finite at all times. This is because, when one follows the motion of a particle 
colliding with the boundary, it may go around the origin indefinitely, and (f>i (which must be defined as a continuous 
variable, without any 2-k jumps) may increase indefinitely. Thus the time average of the time derivative of this 
quantity does not vanish. 
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We now calculate the three terms in the r.h.s. of (3.7). The first one is the average of twice a 
contribution to the Hamiltonian, which is quadratic in the A'^ variables pir (g is diagonal); since 
the average of a quadratic term in the Hamiltonian is 1/(2/3), the first term in the r.h.s. of (3.7) 
is 



1 



m 



N 



N 



The second term reduces to - {{l/2m)Y,f^^Ti{dg'^f /dTi){p, 
replaces {pitp)"^ I2m by l/\^(3g'^'^{Ti)\. The second term becomes 



Finally, since 



1 / ^ 

I3\fr[ tanhTi 



= -(3N 



1 

-pu 



<ro tanh T 



dS 



(3.8) 

Averaging first on pj<^ 
(3.9) 



Je-(^^{dU/dTi)dS2...dSN 



j e-P^dSidS2...dSN 

dn^^\Ti 



dS^ 



the third term can be written as 

Iti<to ^ dri 

Putting together the contributions (3. 8), (3. 9), and (3.11) gives for (3.7) 

_dn(i)(r 



aTo' 



N 1 ro 



tanhr 



+ T- 



dT 



27ra^ sinh r dr 



After an integration by parts, (3.12) becomes 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



3.2 The Thermal Pressure 

The thermal pressure is defined as 



P(^) = - (1^) (3.14) 



dAj 



where F is the free energy. This expression is appropriate for the canonical ensemble, since F 
is related to the canonical partition function Z by PF = —InZ. 



3.2.1 The Thermal Pressure in the Grand Canonical ensemble 

In the following, we shall also need an expression of the thermal pressure appropriate for the 
grand canonical ensemble. It should be remembered that, for a flat OCP in three dimensions, 
the grand canonical partition function must be defined [7] as an ensemble of systems with any 
number N of particles in a fixed volume and with a fixed background charge density —qrif, (using 
an ensemble of neutral systems, i.e. varying together with N would give a divergent grand 
partition function). In two dimensions, /? times the free energy for a neutral flat system [3] 
behaves as [1 — {Pq^ /A)]N \n N as N ^ oo, and therefore the neutral grand canonical partition 
function diverges if jdq^ > 4. This indicates that, in the present case of a 2D OCP on a 



6 



pseudosphere, a similar divergence might occur for an ensemble of neutral systems, and we 
prefer to use an ensemble with a fixed background (which, furthermore, will be seen to be 
exactly solvable at f3q^ = 2). Thus, the grand partition function H and the corresponding 
grand potential O = —{1/(5) InS are functions of /3, A, (, ni,, where ( is the fugacity. The usual 
Legendre transformation from F to Q and from N to changes (3.14) into 



-'"'-(IS) 



(3.15) 



We assume that, even on a pseudosphere, the grand potential is extensive, i.e. of the form 
= Auj{P,(,ni,). Since uj depends on A through rib = ^bl-A-, equation (3.15) becomes 

pW=-a; + n,|^ (3.16) 
drib 

Note the difference with an ordinary fluid, without a background, for which P^^^ = —uj. 
3.2.2 The p(^) - P^*^) Difference 

For a OCP, the thermal pressure is different from the kinetic pressure. In the case of a 2D 
OCP in a flat disk, in the thermodynamic limit, the boundary becomes a straight line and the 
difference was found to be [1] 

p{d) _ p{k) ^ _27rg2nfo / [n^^\x) - nb]xdx (3.17) 

Jo 

where n^^^ {x) is the density at distance x from the boundary. Using the Poisson equation, one 
can write (3.17) in the equivalent form [8] 

_ p{k) ^ gn6[</.,,rface - 0bulk] (3.18) 

where ^buik and ^surface are the electric potential in the bulk and on the disk boundary, respec- 
tively.^ 

Equation (3.18) can be proven as follows. Either in the flat case, or in the case of a pseu- 
dosphere, let us consider a large disk of area A, filled with a 2D OCP. For compressing it in- 
finitesimally, changing the area by dA < 0, at constant P, N, N^, we must provide the reversible 
work SW = -P^^UA. We may achieve that compression in two steps. First, one compresses 
the particles only, leaving the background behind; the corresponding work is SW^^^ = —P^^^dA, 
since P^^^ is the force per unit length exerted on the wall by the particles alone. Then, one 
compresses the background, i.e. brings the charge qrn, dA from a region where the potential is 
^surface i^to the plasma where the potential is (/)(r), spreading it uniformly; the corresponding 
work is 5W^'^^ = qni)dA[(l/A) J 0(r)dS' — (^surface]) where (f>{r) is the potential at distance r from 
the center. Therefore, 

p(^) - pik) = gn6[0surface '^J 0(r) dS] (3.19) 

We expect ^(r) to differ from ^buik only in the neighborhood of the boundary circle. 

In the case of a flat disk, the contribution of this neighborhood to the integral in (3.19) 
is negligible in the thermodynamic limit, 0(r) can be replaced by the constant (phuik, and one 

®In the original papers [1, 8], the derivations of (3.17) and (3.18) have been done in the case of a neutral system. 
However, these derivations can be easily extended to systems carrying a total non vanishing charge. 
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obtains (3.18). On a pseudosphere, (3.19) [with <^(r) instead of (j){r)] is still valid. However, 
now, in the large-disk limit, the integration element dS = 27ra^sinh rdr makes the boundary 
neighborhood dominant, and we rather have 



p{e) _ p{k) 



qub 



(3.20) 



(ro)-e-^° r <j){T)e''dT 
Jo 

After some manipulations, in the thermodynamic limit, (3.20) can be shown to be equivalent to 

pW _ pW = -2Tra^nbq^ / [n^^^ (a) - rih] ae"" da (3.21) 

Jo 

where we have introduced the variable a = tq — t and n^^\a) now denotes the particle density 
at distance aa from the boundary. Indeed, in (3.21), n^^^T) — nb can be expressed in terms 
of 4>{t) through the Poisson equation A^(r) = —2nq[n^^\T) — Ub]- Since the charge density is 
localized at large r, we can use for the Laplacian A ~ a~'^[d'^ /dr'^ + d/dr]. After integrations 
by parts, (3.20) is recovered. 

In conclusion, (3.17) valid for a large flat disc generalizes into (3.21) on a pseudosphere. In 
the limit a — > oo, a ^ 0, aa = x, equation (3.21) does reproduce (3.17). 

3.3 The Mechanical Pressure 

Choquard et al.[l] have also defined a mechanical pressure, in terms of the free energy F, as 

PM = _ (^] 



(3.22) 

' l3,N,nb 

In terms of the grand potential J7, a Legendre transformation now gives 



(3.23) 



If the grand potential is extensive, of the form Q = ALo{P,(^,nb), (3.23) gives 

= -LJ (3.24) 

The difference P^^^-PC') can be obtained by a slight change in the argument of Section 3.2.2. 
Again, we consider a large disk filled with a 2D OOP of area A, and we compress it infinitesimally, 
changing its area by dA < 0, now at constant (3,N,nb, providing the reversible work SW = 
—p{i^)(lj[^ in two steps. Again, first one compresses the particles only, leaving the background 
behind, and the corresponding work is (51^^^^ = —P^'^^dA. Then, one must withdraw the leftover 
background charge qnb dA, bringing it from the surface where the potential is ^Asurface to infinity 
where the potential vanishes. The corresponding work is SW^"^^ = — qubd A (f)suTf3jce- Therefore, 
for a disk on a pseudosphere, — pC^) = qub'psmia.cc-'' 

In the thermodynamic limit, (/isurface 2T:a?q{n — nb) and 

pM _ p{k) ^ 2Tra\'^nb{n - nb) (3.25) 

This difference vanishes for a neutral system (n = nj,). 

The relations (3.21) and (3.25) between the different pressures obtained here by means of 
electrostatic arguments can also be obtained in a more formal way following Choquard et al. [1] , 
using the dilatation method (doing a change of variable r = tqt in the partition function to 
explicitly show the area A dependence) and the BGY equations to replace the two-body density 
terms that appear in the calculations by one-body density terms. 

''This result is identical with the one obtained by Choquard et al.[l] in the case of a flat disk. However, their 
general formula might make difficulties in two dimensions, because the Coulomb potential — ln(r/L) does not 
vanish at infinity and involves an arbitrary constant length L. These difficulties do not arise on a pseudosphere. 
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3.4 The McLxwell Tensor Pressure 



On a pseudosphere, since the area of a large domain is of the same order as the area of the 
neighborhood of the boundary, all the above definitions of the pressure depend on the boundary 
conditions. In previous papers, a definition of a bulk pressure independent of the boundary 
conditions has been looked for. After an erroneous attempt [4], it has been argued [9, 10] that a 
bulk pressure -F^viaxweii could be defined from the Maxwell stress tensor [11] at some point well 
inside the fluid. The result was 




(3.26) 



That same equation of state holds for the 2D OCP on a plane, a sphere, or a pseudosphere. 



4 EXACT RESULTS AT (3q^ = 2 

When the Coulombic coupling constant is Pq"^ = 2, all the thermodynamic properties and 
correlation functions of the two-dimensional one-component plasma can be computed exactly 
in several geometries [3, 12, 13] including the pseudosphere [4]. In reference [4] the density and 
correlation functions in the bulk, on a pseudosphere, were computed. Here we are interested in 
these quantities near the boundary. In reference [4] the calculations were done for a system with 
a — lnsinh((i/2a) interaction and it was shown that this interaction gives the same results as 
the real Coulomb interaction — lntanh(d/2a), as far as the bulk properties are concerned. The 
argument in favor of this equivalence no longer holds for the density and other quantities near 
the boundary; therefore we shall concentrate on the real Coulomb system with a — lntanh((i/2o) 
interaction. This system was briefly considered in the Appendix of reference [4]. For the sake 
of completeness, we revisit here the reduction of the statistical mechanics problem to the study 
of a certain operator. 



4.1 The grand potential 



Working with the set of coordinates (r, i^) on the pseudosphere (the Poincare disk representa- 
tion), the particle i-particle j interaction term in the Hamiltonian can be written as [4] 



v{dij) = — lntanh(c/ij/2a) = — In 



{zi - Zj)/{2a) 



1 - {ziZj/Aa?) 



(4.1) 



where Zj = rje^'^^ and Zj is the complex conjugate of Zj. This interaction (4.1) happens to 
be the Coulomb interaction in a flat disc of radius 2a with ideal conductor walls. Therefore, 
it is possible to use the techniques which have been developed [14, 15] for dealing with ideal 
conductor walls, in the grand canonical ensemble. 

The grand canonical partition function of the OCP at fugacity ( with a flxed background 
density rib, when Pq^ = 2, is 



S = Co 



1 



N=l 



C.{rk)rkdrkdipk 



n 



n 

i<j 



{zi - Zj)/{2a) 



1 - {ziZj/Aa?) 



(4.2) 



where for = 1 the product Wi^j must be replaced by 1. We have defined a position- 
dependent fugacity C,{r) = C,[l — r"^ j [^ia?'^Y™b'^^ -^qC -^yj^ich includes the particle-background 
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interaction (2.10) and only one factor [1 — r^/(4a^)] ^ from the integration measure dS = 
[1 — r^/(4a^)]~^ dr. This should prove to be convenient later. The e*-^ factor is 



e = exp 



ATTUba^ [ Incosh^ — sinh^ Intanh^ 



(4.3) 



which is a constant term coming from the particle-background interaction term (2.10) and 



In Co 



In cosh2 ^ + sinh2 ^ ( sinh^ ^ In tanh^ ^ - 1 



2 V 



(4.4) 



which comes from the background-background interaction (2.11). Notice that for large domains, 
when To — ^ oo, we have 

47rni,a^ 

(4.5) 



oT"0 + l 



and 



InCo ~ — 



(4.6) 



Let us define a set of reduced complex coordinates Ui = {zi/2a) inside the Poincare disk and its 
corresponding images u* = {2a/ Zi) outside the disk. By using Cauchy identity 



det 



1 



Ui — u 



ii,j)e{i,-,Ny 



Uijiui-u* 



(4.7) 



the particle-particle interaction term together with the [1 — (r^/4a^)] ^ other term from the 
integration measure can be cast into the form 



n 

i<j 



{zi - Zj)/(2a) 



2 AT 



j=i 



1 - {ziZj/Aa?) 
The grand canonical partition function then is 

N 



n[l-(rf/4a2)]-i = det — 



UiU 



(4.8) 



Cn 



(4.9) 



We shall now show that this expression can be reduced to an infinite continuous determinant, 

by using a functional integral representation similar to the one which has been developed for 
the two-component Coulomb gas [16]. Let us consider the Gaussian partition function 



Zq = I VipVip exp 



j i){r)M-^{z, ^')V'(r') d\ d\' 



(4.10) 



The fields ifj and are anticommuting Grassmann variables. The Gaussian measure in (4.10) is 
chosen such that its covariance is equal to* 



{il^{Yi)^l}{rj)) = M{zi,Zj 



1 — UiUj 



(4.11) 



* Actually the operator M should be restricted to act only on analytical functions for its inverse M ^ to exist. 
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where (. . .) denotes an average taken with the Gaussian weight of (4.10). By construction we 
have 

Zq = det(M-i) 



Let us now consider the following partition function 

Z = j VipVijjexp J i;{r)M-'^{z,z')i;{r')cfr(fr' + J C(r)Vi(r)V'(r) 
which is equal to 



and then 



where 



Z = det(M-^ + C) 



det[M(M-^ + C)] = det[l + K] 



K{r,r') = M{z,z')ar') 



ar') 



(4.12) 

(4.13) 

(4.14) 
(4.15) 

(4.16) 



1 — uu' 

The results which follow can also be obtained by exchanging the order of the factors M and 
+ ^ in (4.15), i.e. by replacing C('^') by C('') in (4.16), however using the definition (4.16) 
of K is more convenient. Expanding the ratio Z/Zq in powers of ^ we have 



^ = 1 + E 1^ / n d'rkCirk) (V^(ri)V'(ri) ■ ■ ■ V^(riv)V'(rAr)) 
Now, using Wick theorem for anticommuting variables [16] , we find that 

(Vi(ri)V'(ri) • • • Vi(rjv)V(rjv)) = det M{zi, zj) = det ( _^ _ ] 

Comparing equations (4.17) and (4.9) with the help of equation (4.18) we conclude that^ 

E = Co^ = Codet{l + K) 
Zq 



(4.17) 



(4.18) 



(4.19) 



The problem of computing the grand canonical partition function has been reduced to finding 
the eigenvalues of the operator K. The eigenvalue problem for K reads 



c 



12 \ ^Trn^a^ — 1 



ZZ 

4^ 



$(r') r' dr'd^' = A$(r) 



(4.20) 



For A / we notice from equation (4.20) that $(r) = ^(z) is an analytical function of z. Because 
of the circular symmetry it is natural to try $(2;) = ^e{z) = z^ = r^e^^'^ with £ a non-negative 
integer (the functions z^ form a complete basis for the analytical functions). Expanding 



zz' ^AAa'^J 



(4.21) 



1 n=0 

4a2 



Actually, the determinants Zo and Z are divergent quantities, since the eigenvalues of M (restricted to act 
on analytical functions) are easily found to be 47ro^/(£ + 1), with £ any non-negative integer. However, the ratio 
Z/Zo turns out to be finite. 
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and replacing $^(2;) = in equation (4.20) one can show that is actually an eigenfunction 
of K with eigenvalue 

\i = Ana^Ce'^Bto (^ + 1, 4^60^) (4.22) 
with to = ro/(4a^) = tanh^(To/2) and 

Bt^ {i + l, ATTUba'^) = r (1 - t)4™6a2-ii^ (4.23) 
Jo 

the incomplete beta function. So we finally arrive to the result for the grand potential 

00 

= - In S = - In Co - ^ In (1 + Ana^Ce^Bto {£ + 1, ATrnbO^)'^ (4.24) 
e=o 

with e^' and In Co given by equations (4.3) and (4.4). This result is valid for any disk domain of 
radius aro. Later, in Section 4.3, we will derive a more explicit expression of the grand potential 
for large domains tq — ^ 00. 

4.2 The density 

As usual one can compute the density by doing a functional derivative of the grand potential 
with respect to the position-dependent fugacity: 



nW(r)= 1-^ C(r 



5lnE 



(4.25) 



The factor [1 — (r^/4a^)]^ is due to the curvature [4], so that n^^\r) dS is the average number of 
particles in the surface element dS = [1 — (r^/4a^)]~^ dr. Using a Dirac-like notation, one can 
formally write 



In S = Tr ln(l + K) + In Co = J {r\ln{l + C{r)M) \r) dr + In Co 
Then, doing the functional derivative (4.25), one obtains 



(4.26) 



n(^Hr)= 1-^ C(r)(r {l + Ky'M 



r) = 47ra 1-^ C(r-)G(r,r) 



(4.27) 



where we have defined G(r,r') by^'' G = (1 + K) ^M/^Ajra). More explicitly, G is the solution 
of{l + K)G = M/(47ra), that is 



Gir,r') + Ce^ I G(r",r') 



„"2 



4a2 



47rn6a^ — 1 



dr" = 



1 - 



zz 
4^ 



47ra 



1 - 



zz 
4^ 



and the density is given by 

nW(r) = 47raCe^ 



2 \ 47rn(,a^+l 



G(r,r) 



"'the factor 4na is there just to keep the same notations as in reference [4]. 



(4.28) 



(4.29) 
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Prom the integral equation (4.28) one can see that G{r, r') is an analytical function of z. Thus 
the solution is of the form 

oo 

G{r,r')=J2ae{r')/ (4.30) 



and equation (4.28) yields 



1 f ZZ' \ X 

'^^'"''''^ " i^^^li^j l + 47ra2Ce^Sto(^ + l,47rn6a2) ^^'^^^ 
Then the density is given by 

After some calculation (see the Appendix), it can be shown that, in the limit a — > 00, the result 
for the flat disk in the canonical ensemble [17] 

^ - exp(- w^) V (4 33) 

is recovered, up to a correction due to the non-equivalence of ensembles in finite systems. In 
(4.33), 7 is the incomplete gamma function 

j{£ + l,x)= r t^e-^dt (4.34) 
Jo 

In that flat-disk case, in the thermodynamic limit (half-space), n^^\ro) = ncontact n(,ln 2. 
4.3 Lcirge domains 

We are now interested in large domains tq — 00. In this thermodynamic limit we will show that 
the sums in equations (4.24) and (4.32) can be replaced by integrals. For pedagogical reasons 
we will first consider the case ATmbO^ = 1 in which the calculations are simpler, and afterwards 
deal with the general case. 

4.3.1 The case ATrnbO^ = 1 

In this case the incomplete beta function that appears in equations (4.24) and (4.32) simply is 

B (f^lD- - [tanh2(ro/2)]^+^ 

B^,{e + 1, 1) _ ^ _ (4.35) 

When To —> 00 we have 

(4.36, 

Then the sum appearing in the grand potential (4.24) takes the form 



^ / Ce exp(-4(£ + l)e-"o) ' 



(4.37) 
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where we have used the asymptotic expression (4.5) for e*-^. This sum can be seen as a Riemann 
sum for the variable x = 4{l + l)e~'^^ . Indeed, for large values of tq, the variable x varies in 
small steps dx = 4e~^". The sum (4.37) then converges, when tq — > oo, to the integral 



/ 

^0 



ni, X I 4e~'^o 



(4.38) 



This expression together with equation (4.6) for the constant In Co gives the grand potential in 
the thermodynamic limit tq — ^ oo 



(3n — 



Jo 



ln( 1 + ^^ 1 dx-1 
\ rib x I 



(4.39) 



We notice that the grand potential is extensive as expected. The area of the system being 
A = AttoP' sinh^(To/2) 2± 7ra^e'^°, we find that the grand potential per unit area uj = n/ A is given 

by 



-m 



°^ln( l + dx-l 

\ m X I 



(4.40) 



Similar calculations lead from equation (4.32) to the density n^^\a) near the boundary as a 
function of the distance from that boundary aa = a(ro — r). 



Jo 



1 + 



(,ee ^ 

Ub X 



■ dx 



(4.41) 



After the change of variable xe^ x, this can be written as 

7^(1) (cr) foo xe~^ dx 



rib 



f 

Jo 



xe 



(4.42) 



(Ce/nb) 



+ e 



The average density n = N/A can be obtained integrating the density profile (4.42) or by using 
the thermodynamic relation N = — f3C,{d^ / dQ . We find 



n 
rib 



e ^ dx 



iCe/rib) 



+ e 



(4.43) 



4.3.2 The general case 

With the case 47rn{,a^ = 1 we have illustrated the general procedure for computing the thermo- 
dynamic limit. Now we proceed to compute it in the more general case where 47rn(,a^ has any 
positive value. To simplify the notations let us define a = Airnba^- The main difficulty is to find 
a suitable asymptotic expression of the incomplete beta function 

Bt,{£+l,a) = f\l-tY-HUt (4.44) 

when ^ 1 which is valid for large i. As we have noticed in the previous section the main 
contribution from the sum in £ that appears in the grand potential comes from large values 
of i which are of order e^^ . For these values of I the integrand in the definition of the beta 
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function (1 — t)°'~^t^ is very peaked around t = to and decays very fast when t — > 0. So the main 
contribution to the incomplete beta function comes from values of t near • It is then natural 
to do the change of variable in the integral t = tg — v where with the new variable v the integral 
is mainly dominated by small values of v. Then we have 



Bt,ie + l,a)= /*V - to + ^;)°-ie^i°(*o-'') dv 
Jo 



(4.45) 



Replacing to by its asymptotic value to ^ 1 — "^^ and taking into account that v is small (of 
order e~'^°), we find, at first order in e~'^°, 

Bt,ii + l,a)r^^ria,x) (4.46) 

where we have introduced once more the variable x = 4£e~^" (at first order in e~^" it is the same 
variable x = A{£ + l)e~'^° introduced in the case a = 1) and 

POO 

r{a,x)= y'^-'^e-^dy (4.47) 

Jx 

is an incomplete gamma function. With this result and equation (4.5) the term e'~''Btf^{£+ 1)Ck) 
in the expressions (4.24) and (4.32) appears as a function of the continuous variable x = 4£e~'^ 

e^Si„(^ + l,a)~e"^^ (4.48) 

With this result we can replace the sums for £ in equations (4.24) and (4.32) by integrals over 
the variable x and we find the following expressions for the grand potential per unit area 



1 f 

and the density 



2^^47rn(,a2 ^(47rn6a^ x) 



1 + AnaXe . ^ 



dx } (4.49) 



n(i)((7) = ^e4^"''«'e(4'^"*'«'+i)^ H ' " \ ^ (4.50) 

l + 4vra^Ce^-^°^^^^7^°r^ 

In particular the contact value of the density, that is when cr = 0, is 

ncontact = n^'^ (0) = Ce^™^'^^ / ^ (4.51) 

After some calculation (see the Appendix), it can be shown that, in the limit a — > oo, the result 
for the flat disk in the thermodynamic limit ^contact = n^ln 2 is again recovered. 

An alternative expression for the density which we will also use is obtained by doing the 
change of variable xe*^ — ^ x and introducing again a = AnmijO?' 



n 



rib Jo X s. 



(4.52) 



(Ce"M) 



+ aT{a,xe'~'') 
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Prom this expression it can be seen that in the bulk, when o" — ^ oo and e"*^ —>■ 0, the density is 
equal to the background density, n^^\a) n^. The system is neutral in the bulk. The excess 
charge, which is controlled by the fugacity concentrates as usual on the boundary. 

The average total number of particles N and the average density n = N/A can be computed 
either by using the thermodynamic relation 

N = -PC^ (4.53) 

or by integrating the density profile (4.50) 



N = [ nW((T)d5 = 7raVo H n^^\a) e'" da (4.54) 

Jt<tq Jo 

The two methods yield the same result, as expected, 

n=^= Ce^-'"^^ r . ^(4v^n5a^.x)dx 

A Jo x^'^^ca^ +47ra2Ce^""''" r(47rn6a2,x) 

The ratio of the average density and the background density can be put in the form 

^ = r ^i^l^ (4.56) 

Jo ^ , ^Y(a x) 

As seen on equations (4.52) and (4.56) the density profile n^^^ (a) and the average density n are 
functions of the parameter g = (^e^™''" /n^. Different values of this parameter g give different 
density profiles and mean densities. Depending on the value of g the system can be globally 
positive, negative or neutral. Prom equation (4.56) it can be seen that the average density is a 
monotonous increasing function of the fugacity, as it should be. Therefore there is one unique 
value of the fugacity for which the system is globally neutral. For the case Anriba'^ = 1, we have 
determined numerically the value of g needed for the system to be neutral, n = n;,. This value 
is g = (e/rib = 1.80237. 

It may be noted that, in the case of a flat disk in the grand canonical ensemble, the 2D 
OCP remains essentially neutral (the modulus of its total charge cannot exceed one elementary 
charge q), whatever the fugacity ( might be [18,19]; this is because the Coulomb interaction 
— ln(r/L) becomes infinite at infinity and bringing an excess charge from a reservoir at infinity 
to the system already carrying a net charge would cost an infinite energy. On the contrary, in the 
present case of a 2D OCP on a pseudosphere, the Coulomb interaction (2.8) has an exponential 
decay at large distances, and varying the fugacity does change the total charge of the disk. 

Pigure 1 shows several plots of the density n^^^ (a) as a function of the distance a from the 
boundary (in units of a), for different values of g, in the case a = AirnbO^ = 1. It is interesting 
to notice that for g < 1 the density is always an increasing function of a. Par away from the 
boundary, the density approaches the background density n(, from below. On the other hand 
when g > 1, but not too large, the density profile shows an oscillation: n^^\a) is no longer a 
monotonous function of a. Par away from the boundary, a — oo, the density now approaches 
the background density from above. Pinally, when g is large enough, the density profile is again 
monotonous, now a decreasing function of a. 

The change of behavior as a ^ oo can actually be shown analytically. Let us define u = e' 
Prom equation (4.52) we have 



o — IT 



Ub J Jo f{xuY ^. ,^ ^ 

h al (a, xu) 

9 



1 



dx (4.57) 
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The first term in the integral is always positive. The second term, 1 — (e^"/^), in the limit 
cr ^ oo (u ^ 0) is 1 — (1/(7). If 5 < 1 it is negative, then dn^-^^ /du is negative and n^^\a) is 
then an increasing function of a when a —> 00 as it was noticed in the last paragraph. 

Also, in this case a = 1, when ^ — *■ 00 the density profile (4.42) can be computed explicitly 



1 



rib (1 - e- 
It is clearly a monotonous decreasing function of a. 



(4.58) 



4.4 Relations between the different pressures 

Prom the explicit expressions (4.49) and (4.52) for the grand potential and the density profile, we 
can check the relations between the different pressures obtained in Section 3. The mechanical 
pressure simply is = —uj and it is given by equation (4.49). This expression can be 

transformed by doing an integration by parts in the integral giving 



47ra2 



47ra2xCe^'^"''«'-^ 

00 ax 



r(47rnfea^, re) 



2\2 



1 + AwaXe 



2^^47rn.a^ r(47rnba , x) 



dx + {4irni,a ) 



r.47rn(,a^ 



(4.59) 



By the replacement 



d_ 
dx 



T{47rnba'^,x) 



X 



— Aimba 



2r(47rn6a^, x) 



X 



47rni,a^+l 



(4.60) 



in equation (4.59), one recognizes the expressions (4.51) and (4.55) for the contact density and 
the average density, thus giving 



/3pM = n(^)(0) + ATrnbO^in - n^) 



(4.61) 



which is precisely, when /3(/^ = 2, the relation (3.25) between the mechanical pressure P("^) and 
the kinetic pressure pC^) = (l//?)n(i)(0) obtained in Section 3. 
The thermal pressure is 



piO) = -u,iC,nb)+nb[^^j^ 



(4.62) 



The last term in this equation is given by 



Prib 



dub Atto? 
Making the replacement 



2a' 



-f 

Jo 



a- 



1 + 



ATra%e°'T{a,x) da 



e'^T(a,x) 



dx 



(4.63) 



a 



d_ 
da 



e"r(a,x) 



x^ 



T(a,x) d 



r(a, x) 



X" 



da 

in equation (4.63), one recognizes in the first term the average density n, thus obtaining 

Qnb-^ — = a{2nb — n) — al 

OUb 



(4.64) 



(4.65) 
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where 



Jo 



d 



1 + 



47ra%e°'T{a,x) da 



r(a, x) 



dx 



(4.66) 



So the thermal pressure is given by 

pp{e) =„(i)(o) + an6-a/ 



(4.67) 



On the other hand the integral appearing in the general relation (3.21) between the thermal 
pressure and the kinetic pressure 



roo 

J= / {n^^>{a)-nb)e-''ada 
Jo 



can be split into two parts 
with 



J = -Ub + I' 



(4.68) 



(4.69) 



n^^\a)ae-'' da (4.70) 
Using the actual integral representation for the density profile given by equation (4.50) yields 



re 

Jo 



, 47raVe"r(a,x) Uo 
1 H ^ ^ ' 



(4.71) 



X'^ 



The integral over a can be cast in the form 



d_ 
da 



f 

Jo 



(4.72) 



By doing the change of variable y = xe'^ one immediately recognizes the integral representation 
of the incomplete gamma function. The above expression is then equal to 



d_ 
da 



Via, x) 



x^ 



Thus we have proven that I' = I and finally we have the relation 

roo 

/3(p(e) _ p(fc) ) = -A-KTiba^ / (n^^^ {a) - rib) e-^ada 

Jo 

which is relation (3.21) in the solvable case (3q^ = 2. 



(4.73) 



(4.74) 



5 CONCLUSION 

In a flat space, the neighborhood of the boundary of a large domain has a volume which is a 
negligible fraction of the whole volume. This is why, for the statistical mechanics of ordinary 
fluids, usually there is a thermodynamic limit: when the volume becomes infinite, quantities such 
as the free energy per unit volume or the pressure have a unique limit, independent of the domain 
shape and of the boundary conditions. However, even in a fiat space, the one-component plasma 
is special. For the OOP, it is possible to define several non-equivalent pressures, some of which, 
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for instance the kinetic pressure, obviously are surface-dependent even in the infinite-system 
hmit. 

Even for ordinary fluids, statistical mechanics on a pseudosphere is expected to have special 
features, which are essentially related to the property that, for a large domain, the area of the 
neighborhood of the boundary is of the same order of magnitude as the whole area. Although 
some bulk properties, such as correlation functions far away from the boundary, will exist, 
extensive quantities such as the free energy or the grand potential are strongly dependent on the 
boundary neighborhood and surface effects. For instance, in the large-domain limit, no unique 
limit is expected for the free energy per unit area F/A or the pressure —{dF/dA)p,N- 

In the present paper, we have studied the 2D OCP on a pseudosphere, for which surface effects 
are expected to be important for both reasons: because we are dealing with a one-component 
plasma and because the space is a pseudosphere. Therefore, although the correlation functions 
far away from the boundary have unique thermodynamic limits [4], many other properties are 
expected to depend on the domain shape and on the boundary conditions. This is why we have 
considered a special well-defined geometry: the domain is a disk bounded by a plain hard wall, 
and we have studied the corresponding large-disk limit. Our results have been derived only for 
that geometry. 

We have been especially interested by different pressures which can be defined for this system. 
It has been shown that the virial pressure P^'"^ (defined through the virial theorem) and the 
kinetic pressure P^^^ (the force per unit length that the particles alone exert on the wall) are 
equal to each other. We have also considered the thermal pressure P^^\ the definition of which 
includes contributions from the background. It should be noted that this thermal pressure is also 
dependent on surface effects, since it is defined by (3.14) and (3.15) in terms of the free energy 
or the grand potential, and the corresponding partition functions include relevant contributions 
from the surface region. The thermal pressure is not equal to the previous ones. We have also 
considered the so-called mechanical pressure P^™-) which differs from the kinetic one only for 
charged systems. General relations among these different pressures have been established. 

One of the referees of the present paper has asked which one of these different pressures 
is the "right" one, i.e. which one would be measured by a barometer. The answer, based on 
the previous paragraph, is that it depends on which kind of "barometer" is used. For instance, 
the measured pressure would not be the same if the barometer, placed on the wall, measures 
only the force exerted on it by the particles alone, or if it also feels the force exerted by the 
background. 

When /3g^ = 2, the model is exactly solvable, in the grand canonical ensemble. Explicit 
expressions have been obtained for the grand potential, the density profile, and the pressures. 
The general relations between the different pressures have been checked. 

A bulk pressure, independent of the surface effects, can be defined from the Maxwell stress 
tensor. It is not astonishing that this bulk pressure is different from the previous ones, all of 
which depend on surface effects. 

APPENDIX: THE FLAT LIMIT 

In this Appendix we study the flat limit a ^ oo of the expressions found for the density in 
section 4. We shall study the limit a — >■ oo for a finite system and then take the thermodynamic 
limit and compare to the result of taking first the thermodynamic limit and then the fiat limit 
a ^ oo. Since for a large system on the pseudosphere boundary effects are of the same order as 
bulk effects it is not clear a priori whether computing these two limits in different order would 
give the same results. We shall show that, indeed, the same results are obtained. 
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For a finite disk of radius d = utq, we have in the flat hmit a — ^ oo, d ~ rg. In equation (4.32), 
in the hmit a — ^ oo, the term e*-" given by (4.3) becomes 



/ 2 \ 

i e". (A.1) 



where A'^;, = vrnferg is the number of particles in the background in the flat limit. Since for large 
^0 = f'o / {'^0''^) is small, the incomplete beta function in equation (4.32) is 



ii+l,a)= r e^"--'^ dt ~ r dt ~ (A.2) 

Jo Jo 



Expanding (1 — (r^/4a^))'^'^"'''"^ ~ exp(— vrrifer^) in equation (4.32) we finally find the density as 
a function of the distance r from the center 



(vrnftr 

a^-^^-^'-e-^KWC) + 7(^ + 1, N,) 



n(i) (r) = n,e-^-»^' V j, "^-^^^ (A.3) 



When a — oo the terms for l> N^jin the sum vanish because ^ oo. Then 

nW(r) = n,e-w^ -J^^^ + An«(r) (A.4) 

The first term is the density for a flat OCP in the canonical ensemble with a background with 
E{Nf,) elementary charges {E{Nf,) is the integer part of A^^). The second term is a correction 
due to the inequivalence of the ensembles for finite systems and it depends on whether Ni, is an 
integer or not. If is not an integer 



7(i?(A(.) + l,Ab) 

and if Ni, is an integer 



An^^^(r) = nfc ^ (A.5) 



In any case in the thermodynamic limit Tq oo, Nj, ^ oo, this term An^^\r) vanishes giving 
the known results for the OCP in a flat space in the canonical ensemble [3,17]. Integrating 
the profile density (A.4) one finds the average number of particles. For a finite system it is 
interesting to notice that the average total number of particles A" is 

N = E{Nh) + 1 (A.7) 

for A{, not an integer and 

N = Nb + — (A.8) 

for A^E, an integer. In both cases the departure from the neutral case N = Nf, is at most of one 
elementary charge as it was noticed before [18, 19]. 
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Let us now consider the other order of the Umits. We start with the expression (4.51) for 
the contact density in the thermodynamic hmit in the pscudosphcrc and show that in the hmit 
o ^ oo the value of the contact density reduces to the known expression for a neutral OCP in a 
fiat space at a hard wall [17]. We also show that in that limit the average density is independent 
of the fugacity and equal to the background density n = rib- 

Equation (4.51) can be rewritten as 



contact 



rib Jo ^x'^e-"' + ar{a,x) 



(A.9) 



For large a, the numerator of the integrand in (A.9) has a sharp peak at x = a and can be 
expanded as 



(A.IO) 



In the denominator, using the large a expansion of the incomplete gamma function [20], and 
neglecting 1 with respect to a, we obtain 



1 — erf 



x — a + 1 



(A.ll) 



where ^ 

erf(t) = ^ / e-"' du (A. 12) 

■s/tt Jo 

is the error function. Using (A.IO) and (A.ll) in (A.9) gives 



, 2 



"^^^^^^ ~ r ""^i . — XT (A.13) 

ub Jo f{^r + ^[l-evt{^)] 

For X > a, the first term in the denominator goes to infinity for large a and the integrand goes 
to zero. On the other hand, when x < a, this same first term goes to zero, thus, after the change 
of variable t = {x — a)/ \/2a. 



ncontact ^2r0 ^-^^ dt ^^^^^ 

VTT J-^/W^ 1 _ erf + ^ ) 



Finally, as a — ^ oo. 



Q d eri{t) 

'^contact 



rib J-oo 1 - erf(t) 

This is the known value [17] for the contact density at a hard plain wall for a neutral OCP. 
Following the same lines, equation (4.56) for the average density becomes in the limit a — oo 

n_ /2^o [1 - erf (i)] ^ ^ 



rib yaJ-y/^ l-erf(t) 

The average density is equal to the background density and it is independent of the fugacity. 
Whatever value the fugacity has, the system cannot be charged in the flat case in the thermo- 
dynamic limit. 
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List of figure captions 



Figure 1: The density profile n^^\a) (in units of 725) as a function of the distance from the 
boundary a (in units of a) for different values of the parameter g = Qe/nb'm. the case Aimija^ = 1. 
Prom bottom to top, in full line g = 0.5, 1.5, 2.5, 5.0, 10.0 and in dashed line g = 1 (change of 
behavior between monotonous increasing profile and oscillating profile), g = 1.80237 (globally 
neutral system) and 5 — 00. 
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